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We cross-correlate the new 3 year Wilkinson Microwave Anistropy Probe (WMAP) cosmic mi- 
crowave background (CMB) data with the NRAO VLA Sky Survey (NVSS) radio galaxy data, and 
find further evidence of late integrated Sachs- Wolfe (ISW) effect ta king place at late times in cosmic 
history. Our detection makes use of a novel statistical method 1491 |5(]I| based on a new construction 
of spherical wavelets, called needlets. The null hypothesis (no ISW) is excluded at more than 99.7% 
confidence. When we compare the measured cross-correlation with the theoretical predictions of 
standard, flat cosmological models with a generalized dark energy component parameterized by its 
density, Qde, equation of state w and speed of sound c^, we find 0.3 < Qoe < 0.8 at 95% c.l., 
independently of and w. If dark energy is assumed to be a cosmological constant {w = —1), the 
bound on density shrinks to 0.41 < fioE < 0.79. Models without dark energy are excluded at more 
than 4(T. The bounds on w depend rather strongly on the assumed value of c^. We find that models 
with more negative equation of state (such as phantom models) are a worse fit to the data in the 
case Cs ~ 1 than in the case — 0. 

PACS numbers: 98.80.-k, 98.70.Vc, 95.90.+V, 95.75.Pq 



I. INTRODUCTION 

The most outstanding problem in modern cosmology is 
understanding the mechanism that led to a recent epoch 
of accelerated expansion of the universe. The evidence 
that we live in an accelerating universe is now compelling. 
The luminosity distance at high redshift (z ~ 1) mea- 
sured from distant type la supernovae is consistent with 
a negative deceleration parameter (qo < at ~ 3cr) and 
shows strong evidence of a recent transition from decel- 
eration to acceleration Q, |^ Q . The amount of clustered 
matter in the universe, as detected from its gravitational 
signature through a variety of large scale probes (red- 
shift surveys, clusters of galaxies, etc.) cannot be more 
than ~ 1/3 of the total content of the universe P|. Ob- 
servations of the cosmic microwave background (CMB) 
anisotropy have constrained the value of cosmological pa- 
rameters with outstanding precision. The recent WMAP 
data ([SIEI3) have shown that the total density of the 
universe is very close to its critical value. Taken together, 
these results are a strong indication in favor of a non-null 
cosmological term, which would at the same time explain 
the accelerated expansion of the universe and provide the 
remaining 2/3 of its critical density. 

The precise nature of the cosmological term which 
drives the accelerated expansion, however, remains mys- 
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terious. The favoured working hypothesis is to consider 
a dynamical, almost homogeneous component (termed 
dark energy) with negative pressure (or, equivalently, re- 
pulsive gravity) and an equation of state w = p/ p < — 1/3 
[S^ 9, 10, llj. Such a framework helps alleviating a num- 
ber of fundamental problems arising when a constant cos- 
mological term is interpreted as the energy density of the 
vacuum [T^ . 

One key indication of an accelerated phase in cosmic 
history is the signature from the integrated Sachs- Wolfe 
(ISW) effect [13 in the CMB angular power spectrum. 
This is directly related to variations in the gravitational 
potential: in particular, it traces the epoch of transi- 
tion from a matter-dominated universe to one dominated 
by dark energy. This effect (which is usually called late 
ISW, as opposed to a early ISW generated during the 
radiation-matter transition), shows up as a contribution 
in the low multipole region of the CMB spectrum. A 
detection of a late ISW signal in a flat universe is, in it- 
self, a direct evidence of dark energy. Furthermore, the 
details of the ISW contribution depend on the physics of 
dark energy, and are therefore a powerful tool to better 
understand its nature. Unfortunately, the low multipole 
region of the angular power spectrum is also the most af- 
fected by cosmic variance, making the extraction of the 
ISW signal a difficult task. 

A useful way to separate the ISW contribution from 
the total signal is to cross-correlate the CMB anisotropy 
pattern (imprinted during the recombination epoch at 
z ~ 1100) with tracers of the large scale structure (LSS) 
in the local universe 0|. Detailed predictions of the 
ability to reconstruct the ISW using this technique were 
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obtained by a number of authors ITlll7lll9l Is^. 

This kind of analysis has been performed several times 
during the past few years, using different CMB data sets 
and various tracers of clustering. The first detection of 
the ISW 20, 21j was obtained by combining the WMAP 
1st year CMB data with the hard X-ray background 
observed by the High Energy Astronomy Observatory- 1 
satellite (HEAO-1 22]) and with the radio galaxies of the 
NRAO VLA Sky Survey (NVSS HI). The positive cor- 
relation with NVSS was later confirmed by the WMAP 
team 24]. Other large scale structure tracers that led 
to similar positive results were the APM gala xy survey 
m, the Sloan Digital Sky Survey (SDSslP)) and the 
near infrared 2 Micron All Sky Survey eXtendend Source 
Catalog (2MASS XSC [l^) [21121113,111 El El- 

A somewhat different strategy to attack the problem 
was recently adopted by other authors, who attempted 
to seek the ISW signal in spaces other than the pixel 
space of the maps or the harmonic space of the angular 
power spectrum 0, l^l] . This approach relies on spher- 
ical wavelets as a tool to exploit the spatial localization 
of ISW (at large angular scales) in order to get a more 
significant detection of the effect. 

The purpose of the present paper is twofold. On one 
side, we want to perform a further analysis of the CMB- 
LSS cross-correlation, in order to obtain an independent 
check on previous results. We combine the recent 3rd 
year release of WMAP CMB sky maps with the radio 
galaxy NVSS catalogue, and carry out our investigation 
in wavelet space. We make use of a new type of spherical 
wavelets, the so-called needlets 48^ .49^ ^] which have 
a series of advantages over previously used wavelets, as 
will be described in detail later. This then represents at 
the same time a check on previous results [M 111] and 
a significant improvement of the statistical and techni- 
cal aspects of the problem. On the other side, we fol- 
low a rather general approach to dark energy modeliza- 
tion, as first proposed in '36] . Within this framework the 
phenomenology of dark energy is characterized by three 
physical parameters: its overall density r^DE, its equation 
of state w, and the sound speed cj.. This parameteriza- 
tion has the advantage of being model independent, al- 
lowing one to encompass a rather broad set of fundamen- 
tal models, and of giving a more realistic description of 
the dark energy fluid, for example accounting for its clus- 
tering properties, a feature that was shown to have quite 
a strong effect on theoretical predictions IH. As shown 
in 0, [sl m the ISW signature can in principle be able 
to set constraints on the parameters of this generalized 
dark energy scenario: however, previous analyses of the 
ISW from CMB-LSS cross-correlation made a number 
of unrealistic simplyfing assumptions on the dark energy 
component and were only able to either found confirma- 
tions for its existence by constraining its density, or to 
set limits on its equation of state under restrictive hy- 
potheses on its clustering properties (one notable excep- 
tion being the analysis performed in which applied 
a parameterization similar to ours to make a likelihood 



analysis of the cross-correlation data points estimated in 
|40|). Our approach is more ambitious, as we attempt a 
more realistic description of dark energy and derive con- 
straints on the combined set of three above mentioned 
parameters. 

The paper is organized as follows. In Section ^1 we de- 
scribe the dataset we use to perform our analysis. In Sec- 
tion mil we outline the theoretical predictions of the ex- 
pected ISW signal and discuss the statistical techniques 
we apply to extract it from the data. In Section Hvl we 
present our results and the derived constraints on dark 
energy. Finally, in Section^we discuss our main findings 
and conclusions. 



II. DATA 

We trace the local distribution of matter in the uni- 
verse by using the NVSS radio galaxy catalogue [23^ . 
This dataset contains roughly 1.8 x 10^ point sources ob- 
served at 1.4 GHz. The flux limit of the catalogue is at 
~ 2.5 mJy, resulting in a completeness of about 50%. 
The survey covers about 80% of the sky, a.t 6 > —40°. 
We construct a point source map (after removal of about 
3 X 10^ resolved sources) using the iVgide = 64 HEALPix 
pixelization Such a map has 49 152 pixels of about 
1 degree side, and guarantees a good sampling of source 
counts in each pixel. We conservatively exclude from 
the map all sources with S > —37° since the coverage 
becomes very poor when approaching that value of decli- 
nation. The final map we use has roughly 35 sources per 
pixels on average. It was pointed out |4j] that there is a 
declination dependence of the mean source density in the 
catalogue, since the survey had different integration time 
in some well-defined constant-declination bands on the 
sky. As suggested in previous analyses 0, we cor- 
rect for this spurious effect by subtracting the average 
source count in each constant-declination band. 

Our CMB dataset consists of the internal linear com- 
bination (ILC) temperature map from the 3 year release 
of WMAP [41I] . This map is produced by combining 5 
smoothed temperature maps, with weights chosen in such 
a way to produce minimal Galactic foreground contami- 
nation while mantaining the CMB signal. According to 
the WMAP team the ILC map gives a reliable 
estimate of the CMB signal at low multipoles with neg- 
ligible instrumental noise. We believe this is appropriate 
with respect to our goals, since the late ISW effect is ex- 
pected to peak exactly at such large angular scales. As 
an additional caution, we mask out the Galactic plane 
region of the map and bright point sources using one of 
the templates produced by the WMAP team, namely the 
conservative KpO intensity mask 0, ^1 . 

While the original map ILC was produced at a reso- 
lution iVside = 512 in the HEALPix pixelization scheme 
HI] (consisting of 12N^.^^ = 3 145 728 pixels) we degrade 
it to a resolution of A^sidc = 64 to match the resolution of 
the NVSS. This resolution is appropriate for the CMB as 



3 




FIG. 1: The function dN/dz used for the sources in the NVSS 
in our analysis. The dotted curve is the theoretical model 
from 45] , which has a spurious feature due to the breakdown 
of the fit at low z. The continuous line is the fit adopted in 
our analysis. 



which enters in the two-point auto-correlation function 
of the CMB as: 



J2(^LLllBl,crP,icosa) (3) 



An 



where Pi are the Legendre polinomials, a is the angu- 
lar separation between two given points, and the func- 
tion Bt.i models the experimental beam response and 
the pixel window function of the map. 

In an equivalent way, given a projected source count 
map: 



6n 



f]N 

dzb{z)—5{z) 
dz 



(4) 



(where 5 is the underlying matter fluctuation in a given 
direction, b is the bias parameter, and dN/dz was dis- 
cussed previously) we can define the source count auto- 
correlation function: 



well, since we are not interested in the fine-scale details 
of the map. 

A joint mask, including both the KpO mask and the 
NVSS declination limit, is applied to both maps used in 
the analysis. 

There is no redshift information for the individual 
sources in the catalogue. Nonetheless, some knowledge 
of the dN / dz function is needed to connect the observed 
source count fluctuation 5n to the underlying matter fluc- 
tuation 5p (as we will show later). We then use a fit to 
the dN/dz estimated in (4^] and already applied to pre- 
vious analysis of the NVSS catalog j2^. Since the fit in 
[45j breaks down at low redshifts, we have approximated 
it with a Gaussian dN/dZ centered around z ~ 0.9 with 
a width Az ~ 0.8, and normalized in order to give a 
unit integral. The resulting dN/dz used in our analysis 
is shown in Figure^ We have verified that the difference 
at low z has negligible effect on the final results. 



III. TECHNIQUES 

A. Cross-correlation between CMB and LSS 

It is common practice to expand the map of the CMB 
temperature fluctuations into spherical harmonics (Y/m) 
as: 

6T = J2'^LYim{0,cf>) (1) 

Im 

in order to extract the angular power spectrum: 

Cr = {\aL\') (2) 



C^^(a) = {5ni5n2) = 

= T.^^^B%,CrPiico.a). (5) 
/ 

Finally, the cross-correlation between CMB and source 
counts is defined as: 

C™(a) = {ST,6n2) = 

= Y.^^^BT,iB^,CrPi{cosa) (6) 
I 

with the usual definition 

Cr ^ {aLaD (7) 

The theoretical auto and cross-correlation functions 
in a given cosmological model can be calculated by nu- 
merically integrating the Boltzmann equation for photon 
brightness coupled to the other relevant equations, in- 
cluding the linear evolution of matter density perturba- 
tions and the evolution of gravitational potential fluctu- 
ations. We did this by suitably modifing the CMBFast 
code (461] in order to output the needed angular power 
spectra. In particular, we can write the angular cross- 
spectrum in terms of CMBFast temperature and matter 
transfer functions (T; and Ni) as: 

Cf^^47r J fA'{m{k)Ni{k) (8) 

where A'^{k) = k^P{k) /2t:'^ and P(k) is the primordial 
power spectrum of fluctuations. Our CMBFast modifica- 
tion also includes a full treatment of a generalized model 
of dark energy, to be described elsewhere 01 • 
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FIG. 2: Plots of the needlet basis ^pjk for an arbitrary value of 
k and three representative values of j: continuous line j = 7, 
dotted line j = 8, dashed line j — 9. Note that the angle in 
the X axis is measured on the geodesic, not on the tangent 
plane. 



B. Spherical needlets 

In this paper we apply for the first time to cosmologi- 
cal data a new construction of spherical wavelets, the so- 
called needlets. Needlets were introduced into the func- 
tional analysis literature by |48|| ; the investigation of their 
properties from the probabilistic point of view is due to 
|49l |5^ , where their relevance for the statistical analy- 
sis of random fields on the sphere is pointed out for the 
first time. In short, the construction can be described as 



follows. Let Xj = {^jk} 



fe=l,2,. 



denote a pixelization of 



the unit sphere S^; in our case we take the pixelization 
provided by HEALPix • This pixelization divides the 
spherical surface in pixels of the same area, and it is hi- 
erarchical, in the sense that one can build pixclizations 
at higher and higher resolutions, starting from a base 
division of the sphere. The lowest possible resolution 
has 12 pixels; each subdivision increases the number of 
pixels in the map by a factor of 4. The level in the pix- 
elization hierarchy is denoted by the symbol iVgido, with 
^sidc — 1, 2, 4, 8, .... For each level in the pixelization hi- 
erarchy, one has 12 N^^^^ pixels. The pixelization adopted 
sets a limit to the highest multipole number that can be 
reliably extracted from a map. In our case this is given 
by / = 2iV,idc. We then define ;„,ax = [s^'+^] < 2iV,ido 
(with [•] denoting the integer part and s > 1) and choose 
an appropriate value of A^sidc for the j value we want to 
probe. 

Let then ^'(•) denote an infinitely differentiable func- 
tion such that 



5'(u) = for u < — or u > s , 

s 

< *(u) < 1 for i < u < s ; 

s 



(9) 



The numerical implementation of such a function is de- 



scribed elsewhere |50j. Up to a constant factor, we can 
then construct our wavelet basis as: 



where 9 denotes the angle between the two vectors x, ^jk ■ 
Some representative cases for the wavelet basis tjjjk are 
plotted in Figure |5] The wavelets coefficients are then 
derived as 
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(11) 
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here again we have used [•] to denote the integer part of 
a real number. 

In our view, the needlets coefficients f3jk enjoy a num- 
ber of very important properties which are not shared 
by other spherical wavelets. As a first example, we stress 
their localization in the multipole space; the contribution 
from the {aim} is analytically derived to be bounded be- 
tween the multipolcs I = [s-'^^] and I = [s-'^^] ■ In con- 
trast with some of the existing literature, here the depen- 
dence is sharp and explicit on the user-chosen parame- 
ter s. A second important property concerns localization 
properties in pixel space; indeed, it is shown in "4^ that 



\'^Pjkix)\ < 



Cms' 



{l + s^d{x,^,k)r^^-' 



for all X e M = 1,2,3, 



(12) 



where d(x,S,jk) denotes angular distance, M is any posi- 
tive integer and Cm is a constant which does not depend 
either on x or on j. To be more explicit, this result im- 
plies that, outside any neighbourhood of radius e around 
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the direction S^jk {N^{^jk), say), we have 
Cm 



Threshold - 0.1 



for all X i Ne{ijk), M = 1,2,3,... (13) 



that is, the decay of the wavelet transform is faster than 
any polynomial s-' . This clearly establish an excellent lo- 
calization behaviour in pixel space. Note that the con- 
stants Cm do depend on the form of the weight func- 
tion and in particular on the value of the bandwidth 
parameter s; typically a better localization in multipole 
space (i.e., a value of s very close to unity) will entail 
a larger value of Cm-, that is, less concentration in pixel 
space for any fixed j . 

The probabilistic properties of the random coefficients 
/3jfc have been established in [50l |: in that paper, it is 
shown that for any two (sequence of) pixels Cjfe: 'Cjfc' such 
that their angular distance is larger than a positive e, for 
all J, we have 



< 



Cm 



for all M = 1,2,3,. 



(14) 



thus proving wavelets coefficients are asymptotically un- 
correlated as j ^ cxd for any fixed angular distance; equa- 
tion H14|l is clearly a probabilistic counterpart of (|13|l . To 
the best of our knowledge, this is the first example of 
such kind of results for any type of spherical wavelets: 
asymptotic uncorrelation (i.e. independence in the Gaus- 
sian case) simplifies enormously any statistical inference 
procedure. In particular, equation l|14|) is used in to 
derive analytically the asymptotic behaviour of a number 
of procedures based on needlets, including tests on angu- 
lar power spectra or tests for Gaussianity and isotropy. 
In a similar manner, it is possible to show that needlets 
coefficients outside the masked regions are asymptotically 
unaffected by the presence of incomplete sky coverage. 

To provide some numerical evidence of this localization 
property of the needlets coefficients we evaluate, for each 
pixel outside the KpO mask, a normalized squared error 
by taking 
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where (3jk is the wavelet coefficient calculated in presence 
of a sky mask. More explicitly, Djk represents a mea- 
sure of the difference between the needlets coefficients 
evaluated with and without partial sky coverage. In the 
presence of good localization properties in pixel space, 
we should expect this difference to be small outside some 
neighbourhoods of the masked regions. This expectation 
is indeed confirmed by our simulation results, which are 
based upon 100 replicated temperature maps. More pre- 
cisely, we observe that in 64% of pixel outside the masked 
region Djk took values smaller than 0.1; this percentage 
rises to 88% when the threshold value is increased to 0.25 
and to 97% when we take 0.5 as a threshold value. To 
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FIG. 3: Effect of the sky cut on the estimated needlet coeffi- 
cients. Pixels k corresponding to values of the estimator Djk 
above a given threshold are shown in black, and are clearly 
concentrated only very close to the masked regions. For this 
figure, j = 11. From top to bottom, the threshold for Djk 
takes the values 0.1, 0.25, 0.5. 



add some visual evidence as a support to our claims, pix- 
els above the indicated threshold are plotted in Figure |31 
their concentration at the boundary of the masked region 
(i.e. next to the galactic plane, and around point sources) 
confirms the rationale beyond our claims. These issues 
shall be further investigated in an upcoming paper; see 

also ElEi. 



As a final advantage, we stress that our construction is 
explicitly devised for the sphere, and as such it does not 
rely on any sort of tangent plane approximation; we view 
this as an important asset both from the computational 
point of view and in terms of the accuracy of statistical 
results. 

Having extracted the needlets coefficients Pjk from the 
CMB and source count maps, the cross-correlation esti- 
mator in wavelet space, Pj, can be calculated simply as: 
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FIG. 4: The wavelet cross-correlation power spectrum /3j of 
the WMAP and NVSS maps. The points represent the signal 
extracted from the real data, with error bars given by Eq. 
1181 The continues line is the average of the cross-correlation 
power spectra obtained when 1000 simulated CMB fiducial 
data sets are correlated with the real NVSS map: this mea- 
sures the level of correlation expected from casual alignment. 
The shaded area is the la dispersion of the simulated spectra. 
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where iVpix(j) is the number of pixels in the pixelization 
scheme (given by iVpix = 12N^^^^). The theoretical pre- 
diction for (}j can be computed from the expected Cf^ 
as: 



Equation^! We chose the value s = 1.5 in the wavelet 
construction for our analysis. The excess signal peaks at 
value 7 < j < 11, corresponding to angular scales be- 
tween 2° and 10°, as expected from theoretical studies 
01 . The error bars were calculated according to Equa- 
tion^] In order to check that the observed signal was 
not produced by casual alignment of sources in the NVSS 
catalogue with the CMB pattern at decoupling, we pro- 
duced 1000 Monte Carlo simulations of the CMB sky with 
an underlying theoretical fiducial lambda cold dark mat- 
ter (LCDM) model corresponding to the WMAP 3 year 
best fit. The resulting maps were processed through our 
analysis pipeline, and the cross-correlation with the real 
NVSS map was calculated for each simulated data set. 
Figure^lshows the resulting average cross-correlation sig- 
nal (continuous line), which is basically zero on all scales. 
The standard deviation of the simulations is also shown 
in the same Figure (shaded area). These errors, calcu- 
lated through the Monte Carlo procedure, are consistent 
with the analytical estimates of Equation 1181 

The cross-correlation signal extracted from the data is 
significantly higher than the expectation value of the sim- 
ulated data. To quantify the statistical significance of the 
detection, we computed the quantity = — 
(/3]™))VA/3|. We found = 29.8. Assuming is 
distributed as a with 12 degrees of freedom (corre- 
sponding to our 12 data points), we can exclude that the 
(3j were produced under the null hypothesis with 99.7% 
confidence. This result is in agreement with previous 
studies which used the same data 0, 0; Ull ■ We stress 
again the fact that in our analysis the very nature of 
needlets guarantees that the correlation between adja- 
cent data points is very low, even in the presence of sky 
cuts. 
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B. Consequences for dark energy models 



It can be easily shown that the analytical expression for 
the dispersion of the estimated cross-correlation power 
spectrum in needlet space is: 



(2/ + 1) 
167r2 
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which, of course, must be only taken as an approximation 
when dealing with real data. 



IV. RESULTS 

A. ISW detection 

In Figure 0] we show the cross-correlation signal in 
wavelet space extracted from the WMAP and NVSS 
data. The data points shown in the Figure were ob- 
tained following the above described procedure, applying 



We have compared the cross-correlation estimated 
from the data with the theoretical predictions in order 
to set constraints on dark energy models. Following the 
phenomenology described in , ,3(i] , we have modeled the 
dark energy component as a fluid characterized by its 
density parameter, fiDE, its equation of state, w = p/p, 
and the sound speed — Sp/Sp. The latter, needs not 
be the usual adiabatic one, but also accounts for an en- 
tropic sound speed, so that ^ p/p. We have considered 
the above three quantities as the free parameters of our 
analysis. All the other cosmological parameters were held 
fixed at the best fit values estimated using the WMAP 
3 year data: Q = 1, Hq = 7^ km/s/Mpc, f^b = 0.042, 
T = 0.088, A = 0.68 (or ag = 0.74), = 0.951 0. 
When modeling the NVSS catalogue, we adopted a con- 
stant bias parameter b = 1.6, a value suggested in lisl 
and already used in previous analysis, e.g. poll2lLl23 . l34j . 

We have computed the theoretical expectation for the 
cross-correlation by using Eauation ll6l and the output of 
our modified CMBFast code 47]. Our calculation fully 



7 



Cs^ = 

0.0 I — ' — ' — ' — I — ' — ' — ' — I — ' — ^ 
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Hde 

FIG. 5: Constraints at 68%, 95% and 99% confidence level in 
the fioE — w plane. The upper panel was obtained under the 
hypothesis that the dark energy speed of sound is = 0; the 
lower panel was obtained for = 1. 



takes into account the clustering properties of dark en- 
ergy. We restricted our analysis to two different values of 
Cg corresponding to the limiting cases describing a scalar 
field behaviour (c^ = 1) and a matter behaviour (c^ = 0). 

The main results of our analysis are summarized in 
Figure El where we plot the joint constraints on the dark 
energy parameters r^DE and w for the two cases of sound 
speed considered here, and in Figure El where we show 
the separate marginalized likelihoods for JIde and w. 

The first conclusion we can draw from our analysis 
is that the evidence for non zero dark energy density is 
rather robust: we find 0.32 < J^de < 0.78 for — I and 
0.36 < flBE < 0.81 for = 0, both at 95% confidence 
level. A null value of Que is excluded at more than 4cr 
(see Figure |B1 upper panel), independently of c^. When 
we model the dark energy as a cosmological constant (i.e. 
we assume the value w = — 1 for its equation of state), 
the bounds on its density shrinks to 0.41 < Qde < 0.79 
at 95% confidence level. 

On the other hand, the constraints on w are strongly 
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FIG. 6: Marginalized likelihood function for the dark energy 
density Qoe (upper panel) and equation of state w (lower 
panel). In each panel, the continuous curve was obtained 
under the hypothesis that = 1, while the dotted curve is 
for ci = 0. 



influenced by the assumed value of c^, because of the 
different clustering behaviour of dark energy (Figure jSl 
lower panel). In particular, we find that the value of the 
equation of state which corresponds to a cosmological 
constant {w — —1) is well within the 95% c.l. when we 
assume = 0. In this case, we can only put an upper 
bound at 95% c.l.: w < —0.54. When = 1 is assumed, 
we find that phantom models are excluded and that the 
cosmological constant case performs comparatively worse 
than models with larger values of w. Our bounds at 95% 
c.l. are —0.96 < w < —0.16. However, we emphasize 
that, for values of Que ~ 0.7, models with w — —1 are 
a good fit to the data, as it is evident from Figure 
(lower panel). In fact, the LCDM WMAP best fit (with 
Qde — 0.76 and w = —1) has = 9.35, with 12 data 
points. The predicted cross-correlation for some dark 
energy models is shown together with our data points in 
Figure [7\ 
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FIG. 7: The cross-correlation data points estimated in our 
analysis (big dots), with their errors, are shown together with 
the theoretical predictions for a standard CDM model (with- 
out dark energy, dotted curve), the best fit under the hypoth- 
esis that = f (short-dashed curve), the best fit under the 
hypothesis that — (long-dashed curve) and the LCDM 
(w = -1) that best fits the WMAP 3 year data (dot-dashed 
curve) . 



V. CONCLUSIONS 

We have analyzed the WMAP 3 year CMB temper- 
ature data, in conjunction with the NVSS radio galaxy 
survey, and found further evidence of a correlation be- 
tween the CMB fluctuation pattern and the local distri- 
bution of matter, consistent with an ISW effect taking 
place at a late epoch of cosmic evolution. When a flat 
universe is assumed (as suggested by CMB observation) 
the detection of a late ISW signature is a strong evi- 
dence in favour of a dark energy component. Our find- 
ings are based on a new construction of spherical wavelets 
that has a number of advantages with respect to previous 



studies. The presence of a correlation between the CMB 
and the LSS is established with a high level of confidence. 

We have also improved the treatment of the dark en- 
ergy component, introducing a more general parameter- 
ization than those used is similar earlier analysis. Quite 
interestingly, we find that although the case for a non 
zero dark energy contribution to the total density is com- 
pelling and robust, the constraints on w do depend on the 
assumed clustering properties of the dark energy compo- 
nent, namely its sound speed c^. Phantom models, and 
also the ordinary cosmological constant case w — — 1, 
perform worse when a quintessence behaviour = 1 is 
assumed. This is due to the fact that there exist models 
with w ^ —0.4 which predict more correlation at smaller 
angular scales (9 ~ 2°). This is an intriguing result, that 
could imply a ISW effect taking place at redshifts as high 
as z ~ 1, earlier than expected in the cosmological con- 
stant case. A similar preference for larger values of w in 
quintessence models was also found in ^3^]. 

Whether this is an indication of interesting physics 
taking place between the dark energy and dark matter 
components is a subject that requires further investiga- 
tion. Clearly, the observation of ISW is proving quite 
promising as a tool to answer the questions arising from 
the mysterious nature of dark energy. While the CMB 
data have reached a great degree of accuracy on the an- 
gular scales that are more relevant for the detection of 
ISW, deeper redshift surveys and better catalogues can, 
in the future, improve the tracing of the local matter 
distribution, thus allowing to reduce the errors on the 
cross-correlation determination. 
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